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Abstract 

A new category of front propagation problems is proposed in which a spread- 
ing instability evolves through a singular configuration before saturating. We 
examine the nature of this front for the viscous Rayleigh instability of a col- 
umn of one fluid immersed in another, using the marginal stability criterion to 
estimate the front velocity, front width, and the selected wavelength in terms 
of the surface tension and viscosity contrast. Experiments are suggested on 
systems that may display this phenomenon, including droplets elongated in 
extensional flows, capillary bridges, liquid crystal tethers, and viscoelastic flu- 
ids. The related problem of propagation in Rayleigh- like systems that do not 
fission is also considered. 
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Front propagation problems are often divided into two broad classes: the invasion of 
one stable state by another, and the invasion of an unstable state by a stable one [0. Here 
we point out the possibility of a third, qualitatively different type of front propagation: the 
invasion of an unstable state by one that develops discontinuously into its final configuration. 
Although we expect this class to be very broad, we focus on one particular example, the 
Rayleigh instability, wherein a cylindrical body of fluid longer than its circumference is 
unstable to breakup into droplets 0. This is a linear instability || triggered by infinitesimal 
perturbations of sufficiently long wavelength. 

Recent work has shown Rayleigh instabilities in many experimental systems: laser- 
tweezed cylindrical lipid vesicles [|J, liquid crystal tethers embedded in a polymer matrix 
||, shrunken polymer gels ||, and columns of viscoelastic fluid 0. Traditional experiments 
display an instability developing uniformly along the cylinder M. However, experiments 
on the "pearling instability" of laser-tweezed membranes reveal a different scenario — front 
propagation. The peristaltic shape deformation spreads out from the laser spot with a con- 
stant velocity . The emerging explanation |10|-[P2]| for this instability is that it is driven by 
the tension induced in the membrane by the laser trap. Fissioning is ultimately prevented 
by the membrane elasticity. 

Motivated by the phenomenology of the pearling instability, we ask, "Can the Rayleigh 
instability develop as a propagating front?" (as in Fig. |l|). Propagation in this context 
is problematic. As a droplet pinches off, the two tips of the broken neck recede from the 
pinching point; if the retracting neck overtakes the front, propagation will be spoiled. Classic 
experiments on the breakup and relaxation of elongated droplets Jl3| reveal that this com- 
petition depends on the viscosity contrast of the two fluids. This retraction and topological 
rearrangement of interface is absent in overdamped systems known to have propagation, 
such as dendritic growth, viscous fingering, and reaction-diffusion systems. The mathe- 
matical structure of our problem differs from those just mentioned as well, with a global 
conservation law from fluid incompressibility, nonlinearities from the two principal surface 
curvatures, and singular behavior near pinching. The former two features are also common 
to Rayleigh-like systems that, like membrane tethers, do not fission. 

The novel aspects outlined above are the very features that complicate both analytical 
and numerical analysis of our problem. As a first step, we apply the intuition gained from 
front propagation problems in local partial differential equations (PDEs) to the Rayleigh 
problem, with the hope that these estimates will serve as a guide for further investigation. 

With the exception of the one-dimensional nonlinear diffusion equation [III], there are 



no rigorous analytic results for front propagation in PDEs. On the other hand, in many 
cases the marginal stability criterion (MSC) |T3|] correctly predicts the front properties. We 
apply the MSC to the Rayleigh instability to make concrete predictions for front velocity, 
front sharpness, and selected wavelength, and suggest some experimental systems where this 
phenomenon may be seen. 

The three categories of front motion are illustrated by generalized Fisher-Kolmogorov 
(FK) equations [17]], 

where £ — J dx{u 2 x /2 — V(u)), with three types of potentials shown in Fig. 0. Recall the 
standard mechanical analogy for traveling-wave solutions: the ansatz u(x,t) = f(x — vt) 
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reduces (JT]) to Newton's equation for a particle of coordinate u at time z = x — v t in a 
potential V(u) with friction coefficient v. In case (a) there is a unique v that allows the front 
to connect the two locally stable states. When a stable state invades an unstable one (b), the 
particle reaches u = at z = oo for every friction coefficient v > 0. However, for sufficiently 
localized initial conditions, there is a unique front velocity that appears asymptotically, and 
is correctly predicted by the MSC. Finally, (c) displays a situation loosely analogous to the 
new category. A potential of this form, e.g. V(u) = u 2 /2 + w 4 /4, leads to a finite-time 
blow-up of u ||l8fl . Although this model does not describe a topological transition, it is like 
the Rayleigh problem in which the growing instability has a finite-time singularity. 

This mechanical analogy suggests that case (c) would, like (b) , possess a family of possible 
velocities. With this motivation, we examine the Rayleigh problem using the linear MSC, 
which yields predictions directly from the linear growth rate u(q) The front velocity 
v* and the shape of the front's leading edge, u(x, t) ~ exp(u;(g*)t + iq*x), are given by the 
relations |T^[ 

„ Re^* du* dcu* 
v = =— -, Im^r- = 0, v = Re—-, 2 
Reg* oq oq 

where u* = w(q*) and q* = q' + iq" is complex. In the standard MSC treatment, the 
selected wavenumber go of the saturated pattern is different from the selected wavenumber 
q' in the leading edge of the front; go is deduced by working in the rest frame of the front 
and assuming nodes are conserved as they pass through the front and into the saturated 
pattern. We do not expect go to be relevant to the Rayleigh problem since the evolution of 
the interface is not continuous. 

The conservation of fluid in the Rayleigh problem is like that appearing in the invasion of 



the disordered phase of a block copolymer by the lamellar phase pD |, with global conservation 
of the mean composition. In contrast to the FK equation, this leads to a dynamics of the 
form u t + j x = 0, with the flux j = —d x (5£/5u), or 

— — • (3) 

dx 2 5u 

Consider now a thread of fluid of radius R and viscosity i]~ in another fluid of viscosity 
i] + = i]~/a, with interfacial tension E. The growth rate of a volume-preserving axisym- 
metric sinusoidal perturbation with wavenumber q = k/R follows from the Stokes equations 
^ ± V 2 v ± = ^p ± , subject to the boundary conditions of no-slip, no tangential stress, and 
the Laplace law for the jump in normal stress, (cr^ — cr^nj = —TLHrii. Here, Hi is the out- 
ward surface normal, H is the mean curvature 2H = r xx (l + r 2 )~ 3 ^ 2 — r _1 (l + r 2 )" 1 ^ 2 , 
and the stress tensors are af* = ^{SJivf + V jvf ) — p ± 5ij. The growth rate is |21 



cu(k,a) = (T,/Ri] + )A(k,a)(l — k 2 ) . The dynamical factor A(k,a) (too lengthy to quote 
here) accounts for the dissipation, and the factor S(l — k 2 ) is associated with the energy of 
a distortion. 

We gain intuition about A(k, a) from its small k form: 



A 



-k 2 


-§ + a(\ + | 7 - \ log 2) + \a log k 




a — k 2 (a 


-1) 


J + o:(f + § 7 -flog 2) + fa log A; 
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where 7 is Euler's constant. For finite a the asymptotic limit as k — > is A(k, a) ~ 
— \k 2 log k. The logarithm is due to large-scale flows in the outer fluid, while the two powers 
of k are analogous to the two x-derivatives in ([3]). They can be understood intuitively by 
noting that one power comes from the flux form of the dynamics, nd t r 2 + d x J = [|l^], with 



flux J ~ r 2 U, where U is the average axial velocity inside the tube. For small k the Stokes 
equations imply a Poiseuille flow of the inner fluid and thus the second power of k through 
Darcy's law for the average velocity, U ~ —r 2 p x . 

If we consider the case of a = 00 (no outer fluid), we recover the small k limit of 
Rayleigh's stability analysis [Q, which has a non-zero growth rate at k — 0. The inner fluid 
has plug flow since it does not have to entrain the outer fluid as the shape of the boundary 
changes. This in turn modifies the relationship between the average velocity U and the 
pressure so for small k, U xx ~ p x |22|| . 

Fig. H displays the results of the MSC applied to Tomotika's growth rate for a range 
of viscosity contrasts a = ri~/r] + |[23|| . For example, the results for the castor-oil-eugenol 
mixture immersed in silicone fluid (a = 0.0315) used in [16[] are the following: the fastest- 
growing mode fc max = 0.502, k! = q'R = 0.404, k" = q"R = 0.291, and v* = 1.04£/r/+ The 
MSC predicts a relatively sharp front; if we measure the front width in units of 2ir/k' we 
see that it is generically a fraction of a growing bulge, from around 0.25 at a = 1 to about 
0.1 for large and small a. We note also that v* differs markedly from the naive estimate 
u;(/c max , a)/fc max . From these results we estimate the distance (k'/2n)(v */oj{k')) the front 
moves in the time it takes a bulge to pinch off; it is about twice the front width for the range 
of a shown, so the droplets pinch off right behind the front. 

Now we crudely estimate the retraction velocity of the bulge at the end of a broken 
thread to see when we expect to find propagation. Assuming a spherical bulge, we balance 
the tension force with the drag of a liquid sphere of radius R to deduce a velocity 
V = (£/?7 + )(r/i?)(a: + l)/(2a + 3), where r < R is the radius of the neck attaching the bulge 
to the rest of the thread. Note that V does not depend sensitively on the inner viscosity. 
Comparing this expression with Fig. 3, we see that the front will outrun retraction at small 
a but not for large a. Therefore, we expect the MSC to describe front propagation at small 
a if it applies at all. 

These predictions are consistent with the experimental results of |13| , which concern the 
relaxation and breakup of an initially elongated drop in an otherwise quiescent fluid. It 
is found that for large a the bulbs on the ends of the dog-bone shaped droplet retract a 
substantial amount before they pinch off (if they pinch off at all); for lower values of a, 
the retraction is relatively slower and the ends do not move as much before they pinch off. 
As alluded to above, this behavior occurs because the outer viscosity limits the retraction 
velocity and the inner viscosity inhibits pinching ||13|| . The experimental results at low a are 
tantalizingly close to our conception of a propagating Rayleigh instability. We note that the 
mechanism of "end-pinching" in which the end of the extended drop bulges and pinches off 
while the rest of the drop is stationary |13| is consistent with the MSC prediction of a sharp 
front. It remains to be seen if these disturbances spread at a constant velocity. 

Since the applicability of the MSC to our problem is speculative, it is important to 
check our predictions. A numerical computation of the interface evolution in the fissioning 
Rayleigh problem is a challenging task. As a first step to test the internal consistency of 
the MSC we studied a model PDE for a Rayleigh-like instability that approaches pinching, 
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with an energy functional with an elastic contribution like that for membranes. Using a 
Poiseuille flow approximation to relate the fluid flux to the pressure gradient, we obtain 




(4) 

with S = J dS {£ + (1/2)k(2H) 2 } and k the bending modulus. Note the natural progression 
from (|I|) to (|3|) to Q). Although (^) has a non-singular evolution, it has the conservation laws 
and nonlinear structure emphasized in the introduction. This structure is similar to thin- 
film dynamics [^|| in which propagating Ray leigh- Taylor instabilities have been observed. In 
T2fl we found that supported front solutions; Fig. f| shows close agreement between the 



numerical solution of ([|) and the MSC The discrepancy is likely due to the slow approach 
of the velocity to its asymptotic limit [|15[ coupled with computational limitations. Note in 



Fig. |] how as each bulge forms it is pulled away from the thin tether, in accord with our 
argument about tension and drag. 

We now turn to possible experiments. The examples we list here do not necessarily have 
a singular evolution; they illustrate the other novel features of propagation outlined in the 
introduction and raise as well the issue of the mechanism by which pinching is inhibited 
in each case, (i) In addition to droplets in extensional flow, another ideal candidate for 
propagation with pinching is the neutrally buoyant capillary bridge stabilized by an electric 
field |16| . (ii) Mather et al. || are studying a shape instabilities of a thread of liquid crystal 



polymer embedded in a polymer matrix, with a competition between tension and nematic 
elasticity, (iii) Matsuo and Tanaka || have stretched a cylindrical gel between two glass 
plates and observed a peristaltic instability in the presence of an appropriate solvent. In 
addition to the long time scales, these systems have the advantage that thermal fluctuations 
are small. These facts may conspire to allow one to see a propagating Rayleigh instability 
in the cases without a precisely controllable tension, (iv) Jets of viscoelastic fluid, as in 
the work of Renardy, et al. jTJ, may also display propagating Rayleigh instabilities, (v) 
Finally, we expect there are examples of front propagation with singular evolution that are 



not Rayleigh-like. For example, in surface-tension-driven Benard convection |[27|| ) an initially 
flat layer of liquid distorts to expose dry spots when it is heated from below. We suggest that 
these dry spots can spread behind a front that moves at constant velocity, rather analogous 
to those seen in the Rayleigh- Taylor instability |25 . 



We have argued that front propagation in Rayleigh and Rayleigh-like systems is qual- 
itatively different from more familiar examples of propagation. The differences are both 
mathematical and physical: the non-pinching Rayleigh-like systems have global conserva- 
tion laws and a particular nonlinear structure due to the two dimensional nature of the 
interface, whereas the Rayleigh problem in addition has singularities and the new physical 
ingredient of retraction. If these instabilities do indeed propagate, it will be important to 
determine whether the propagation can be described by the MSC calculations presented 
here or if some new picture is needed. There is a clear need for future experimental and 
theoretical investigations. 
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FIGURES 



FIG. 1. Propagation of the Rayleigh instability. 

FIG. 2. Potentials illustrating different classes of front propagation in the mechanical analog 
for traveling-wave states. 

FIG. 3. Length scales and front velocity from the marginal stability criterion. 

FIG. 4. Numerical results. Propagation of the Rayleigh instability with an elastic cutoff (top). 
Front velocity as a function of tension (bottom). Solid line is MSC prediction. Dashed line indicates 
critical tension for onset of instability. 
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